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ABSTRACT 

The COMPTEL instrument on the Compton Gamma Ray Observatory is used to 
measure the locations of gamma-ray bursts through direct imaging of MeV photons. 
In a comprehensive search, we have detected and locahzed 29 bursts observed between 
1991 April 19 and 1995 May 31. The average location accuracy of these events is 1.25° 
(Icj), including a systematic error of ~0.5°, which is verified through comparison with 
Interplanetary Network (IPN) timing annuli. The combination of COMPTEL and 
IPN measurements results in locations for 26 of the bursts with an average "error 
box" area of only ~0.3 deg^ (1<t). We find that the angular distribution of COMPTEL 
burst locations is consistent with large-scale isotropy and that there is no statistically 
significant evidence of small- angle auto-correlations. We conclude that there is no 
compelling evidence for burst repetition since no more than two of the events (or ~7% 
of the 29 bursts) could possibly have come from the same source. We also find that 
there is no significant correlation between the burst locations and either Abell clusters 
of galaxies or radio-quiet quasars. Agreement between individual COMPTEL locations 
and IPN annuli places a lower limit of ~100 AU (95% confidence) on the distance to 
the stronger bursts. 

Subject headings: gamma rays: bursts, observations — methods: statistical, numerical 
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1. Introduction 

Cosmic gamma-ray bursts (GRBs) are among the brightest, most plentiful astrophysical 
phenomena in the observable gamma-ray sky, yet they present unique observational challenges. 
The greatest challenge is in measuring burst locations accurately. Because of this difficulty, 
much of our current understanding has come from large-scale statistical measures, where the 
accuracy of individual burst locations is not as important as the total number of bursts. The 
Burst and Transient Source Experiment (BATSE) has now measured the locations (accurate to 
^^5°, on average) and intensities of more than 1700 GRBs — yielding evidence that the sources 
are highly isotropic on the sky, yet radially inhomogeneous (e.g., Meegan et al. 1992; Briggs et 
al. 1996; Pendleton et al. 1996). These results rule-out the once-favored galactic neutron star 
paradigm (e.g.. Hurley 1986; Higdon & Lingenfelter 1990), but still allow three vastly different 
distance regimes: heliocentric (distance D < 10^ AU), extended Galactic halo [D > 200 kpc) 
and cosmological (D ~ 1 Gpc). The GRB mystery is thus far from being solved, as the range 
of possible distances spans over 12 orders of magnitude and the source objects remain elusive. 
In this paper we use the COMPTEL instrument aboard the Compton Gamma Ray Observatory 
(CGRO) to measure GRB locations with improved accuracy. A smaller sample of more-accurate 
burst locations cannot improve upon large-scale measures, but improved location accuracy is 
a fundamental requirement for important studies of small-scale angular clustering, correlations 
with known objects and searches for counterpart objects at other wavelengths. In addition, the 
large-scale distribution of unique sub-populations could provide important insights. 

Recent investigation of small-scale angular clustering of GRB sources has been hampered 
by the relative inaccuracy of locations measured with BATSE. Several statistical studies have 
yielded conflicting results indicating the presence (Quashnock & Lamb 1993; Wang k, Lingenfelter 
1993) and absence (Meegan et al. 1995; Brainerd et al. 1995; Tegmark et al. 1996) of small-angle 
correlations among BATSE burst locations. The presence of small-scale angular correlations 
would be important as it could indicate repeating GRB sources, or unresolved clusters — both of 
which are difficult to explain in most scenarios involving catastrophic destruction or mergers (e.g., 
Meszaros & Rees 1992; Narayan, Paczyhski & Piran 1992; White 1993; Katz 1993; Bickert & 
Greiner 1993). Earlier studies found no evidence to indicate burst repetition on time scales shorter 
than ~10 yr, but a majority of the locations used in these studies were also relatively inaccurate 
(Schaefer & Cline 1985; Atteia et al. 1987). 

The investigation of cross-correlations between GRBs and catalogs of known objects has also 
been limited by imprecise locations. Most studies have failed to reveal any significant association 
with any Galactic or extragalactic source catalog (Howard et al. 1994; Nemiroff, Marani & Cebral 
1994, Harrison, Webber & McNamara 1995). Recently, however, BATSE 3B burst locations have 
been reported to be weakly correlated (~2-3c7 significance) with galaxy clusters (Kolatt & Piran 
1996; Marani et al. 1997; but see Hurley et al. 1997) and radio-quiet quasars (Schartel, Andernach 
k, Greiner 1997). Such correlations might be expected at some level if GRBs lie at cosmological 
distances characteristic of the galaxy samples, where the brighter, closer bursts are expected to 
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show the strongest signal. However, no significant correlations with either galaxy clusters or 
quasars were found in an earlier analysis of accurate locations of bright GRBs, where one might 
expect a stronger signal (Webber et al. 1995). 

The most crucial need for accurate burst locations is in the identification of counterpart objects 
at other wavelengths — a potential key to solving the GRB mystery. Unfortunately, quiescent 
counterpart emission at any wavelength appears to be very weak, or non-existent (see Vrba 1996). 
The last remaining hope in this field is to detect transient counterpart emission soon after the 
onset in gamma rays — a particularly difficult challenge as it requires quick and accurate GRB 
locations. Most simultaneous and near-simultaneous measurements that have been obtained are 
relatively insensitive, but indicate that if low-energy (i.e., non gamma-ray) counterpart emission 
exists at all, it must be weak and/or short-lived (< 1 day). Accurate and timely GRB localization 
is thus crucial for future counterpart search efforts (see McNamara, Harrison & Williams 1995). 
The recent discovery of weak x-ray/optical emission from the directions of GRB 970228 and GRB 
970508 offers the most convincing evidence obtained thus far that GRBs do produce transient 
emission at other wavelengths (Costa et al. 1997a, 1997b; van Paradijs et al. 1997; Bond 1997). If 
these bursts arc typical, then other counterparts will be detected only through the use of sensitive 
instruments combined with rapidly determined, accurate GRB locations. 

The COMPTEL instrument has the unique capability of directly imaging MeV photon 
sources over a wide (--^1 sr, FWHM) field-of-view (FoV). It is thus able to measure the locations 
and energy spectra of several GRBs per year. Full-sky exposure to bursts is achieved through long 
observations, over several years, wherein the CGRO satellite is regularly re-pointed. In this paper 
we use the MeV imaging capability and effective full-sky coverage of COMPTEL to accumulate a 
comprehensive catalog of GRB locations observed during the first four years of the CGRO mission. 
The catalog constitutes a unique sample of relatively accurate locations that is used to investigate 
the GRB spatial distribution — providing insight into the distance/source object mystery. The 
energy spectra and temporal properties of these events have been discussed elsewhere (Kippen et 
al. 1996a, 1997 and references therein), as have preliminary locations of some of the bursts (see 
Table 1). 

2. Instrumentation and Fundamental Analysis 

2.1. Instrument Characteristics 

COMPTEL is a wide-field, double-scatter Compton telescope capable of imaging gamma-ray 
sources in the energy range 0.75-30 MeV. It consists of two planar layers of scintillation detector 
modules separated by ~1.5 m. The upper (Dl) layer incorporates seven low Z, low density liquid 
(NE213A) detectors (total area 4188 cm^), while the lower (D2) layer is composed of 14 high Z, 
high density Nal(Tl) modules (total area 8620 cm^). A telescope event is defined as a coincident 
set of interactions in a single pair of Dl, D2 modules in anti-coincidence with a signal from any 
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of four plastic charged particle veto domes. For each telescope event, the energy deposits and 
interaction locations within each detector are measured using the signals from photomultiplier 
tubes affixed to the detector modules. The interaction sequence from Dl to D2 is identified by 
a time-of-flight (ToF) coincidence measurement and the absolute time of each event is tagged 
with an accuracy of 125 //s. Rejection of a large fraction of neutron-induced background events is 
realized through the measurement and discrimination of pulse-shape (PSD) in Dl. 

The COMPTEL design favors the ideal interaction wherein an incident photon undergoes a 
single Compton scatter in Dl, followed by complete absorption of the scattered photon in D2. For 
such ideal telescope events, the measurement of energy deposits in Dl and D2 yields the scattered 
electron energy and scattered photon energy, respectively, while the measurement of interaction 
locations within the detectors reveals the direction of the scattered photon. Combined, the energy 
and location information constrains the incident photon direction to an annular ring or event circle 
on the sky. The angular radius (f of an event circle is determined by the kinematics of Compton 
scattering: 



where mgC^ is the electron rest energy and Ei and E2 are the energy deposit measurements in 
Dl and D2, respectively. Ambiguity in the arrival direction of individual photons is removed 
through the statistical superposition of many events with different scatter angles, thereby forming 
a unique source signature. Complicating matters are non-ideal source events, such as those arising 
from multiple Compton scatters or pair production in Dl or incomplete energy absorption in 
the detectors, for which equation (|^) is invalid. Non-ideal events result in a complicated source 
signature that must be accounted for in the imaging analysis. Selection criteria applied to event 
parameters (i.e., Ei, E2, (f, ToF and PSD) during analysis can significantly alter the effective 
energy range, sensitivity and FoV of the instrument. A complete description of COMPTEL and 
its operation and performance was presented by Schonfelder et al. (1993). 



COMPTEL employs two independent operating modes for observing transient phenomena 
such as gamma-ray bursts. In parallel with the main telescope mode, a secondary burst mode uses 
two of the D2 modules to accumulate time-resolved spectra in the energy range 0.1-10 MeV (see 
Hanlon et al. 1994 and references therein). In this paper we are concerned only with the source 
localization capabilities of the telescope mode observations. 

In detecting and localizing gamma-ray bursts, COMPTEL is in most cases a photon-limited 
device, whose sensitivity is governed by low detection efficiency rather than background. Sensitivity 
in this case depends on burst fluence, with the minimum detectable fluence given by 




(1) 



2.2. Burst Detection Capabilities 
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where E is the spectrum-dependent mean photon energy (0.75-30 MeV), ^eflf is the effective 
detection area for ideal events and Ns is the number of ideal source events that are required to 
construct a meaningful image (Winkler et al. 1986). Requiring 20 ideal events from a normally 
incident E~'^ power law burst source (Agg = 31 cm^; E = 2.8 MeV) yields S'min = 4.5 x 10~^ 
erg cm~^. Burst sensitivity decreases as the zenith angle t^c between the telescope pointing axis 
and the source direction increases. This relation effectively limits the FoV for bursts to < 60°. 
It is illustrated in Figure la for several different power law energy spectra. The strong dependence 
of burst sensitivity on flucncc (rather than, e.g., peak flux) results in a significant bias against the 
detection of short duration bursts, since short bursts typically do not result in enough events for 
meaningful imaging (even though the instantaneous flux may be large) . 
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Fig. 1. — The estimated sensitivity of COMPTEL to gamma-ray bursts, (a) The minimum 
detectable fluence Smin (0.75-30 MeV) as a function of burst viewing angle i^c for several different 
power law energy spectra, assuming zero background. For bursts longer than ~25 s the estimates 
in (a) must be multiplied by a background correction factor (b), which also depends on the viewing 
angle, and can change by a factor of ~3 due to orbital variations. 
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The average rate of background events that could be confused with ideal source events ranges 
from ~0.2 Hz at = 50° to ~0.6 Hz at -de = 10°, with orbital variations as large as a factor of 
three. Background is thus significant only in the detection of longer duration, low fluence bursts. 
The statistical detection significance for such bursts is 

<7dct ~ , = , (3) 

^N;. + 2Nb 

where A^g is the number of source counts (related to the burst fluence through eq. and A^b is 
the number of background counts (roughly proportional to burst duration). When the number 
of source counts required for a significant detection above background (-/Vg) exceeds the number 
of counts required for imaging (A'^g ~ 20), detection is background-limited. In this case, the 
sensitivity estimates of Figure la are too low. Correction factors (roughly independent of spectral 
shape) for iTdct = 3 are plotted as a function of burst duration in Figure lb. The detection of 
bursts shorter than ~25 s is effectively always photon-limited regardless of zenith angle. For 
longer bursts, background becomes increasingly important. Fortunately, most long duration bursts 
contain short, intense episodes whose detection is less dependent on background. 

For very intense bursts, or burst intervals, where fluence S > 10~^ erg cm~^, electronics and 
telemetry deadtime limit the ability of COMPTEL to register and record events. This affects the 
detection of short, intense bursts, where signiflcant numbers of events fail to be registered and/or 
recorded. Deadtime thus produces an additional bias against the detection of short duration 
bursts. 



2.3. Burst Imaging Analysis 

The imaging analysis of COMPTEL data is performed in a three-dimensional dataspace 
defined by the quantities (x,'i/')<^)- The coordinates (x,^') define the direction of the vector 
between the Dl, D2 interaction locations and the scatter angle (p is computed from the energy 
deposit measurements using equation ([l|). In this representation, the point source signature, or 
point spread function (PSF), is a partially filled cone with 45° opening angle whose vertex is 
centered on the source position (Figure 2). The complicated diagonal shape of the cone mantle 
is determined by the inherent measurement errors and the physics of non-ideal interactions, that 
are all functions of source position and energy spectrum. Dependence of the PSF shape on source 
position is weakened by factoring-out a geometry term Q (see Schonfelder et al. 1993). This allows 
(to some degree of approximation) the use of a single PSF for imaging sources located anywhere 
in the FoV. The PSF used for all burst source analysis in this paper has been derived from Monte 
Carlo simulation (see Stacy et al. 1996) of the COMPTEL response to a point source of gamma 
rays at -de = 10°, with an assumed E^'^ power law energy spectrum. Various standard selections 
are imposed on the event data in order to reduce systematic effects due to incomplete knowledge 
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of the PSF and to improve the signal-to- noise ratio in the data (see Schonfelder et al. 1993). These 
selections are applied equally to the PSF simulations and the data. 



Fig. 2. — Illustration of the three-dimensional COMPTEL response to a point source located 
at sky coordinates (x', '(/'')• Contours show a cross-section of the response function indicating the 
off-diagonal effects of intrinsic measurement errors and non-ideal interactions. 

The imaging process amounts to searching for the PSF signature in the binned (x, V') 
distribution of telescope events (a mixture of source and background). The bin-size is chosen to be 
on the order of the instrumental resolution (1°, 1° and 2° in x, -0 and respectively), resulting 
in few events per bin. A maximum likelihood search technique capable of handling the sparsely 
populated dataspace is employed to obtain quantitative constraints on the source parameters. 
The maximum likelihood method, described fully by de Boer et al. (1992), attempts to fit a 
point source (PSF) plus a model of the background signature to the 3-d dataspace distribution of 
events assuming Poisson counting statistics. For a trial source position, the likelihood of source 
plus background £(S-|-B) and the likelihood of the null hypothesis of pure background /2(B) are 
independently maximized by allowing PSF and background scale factors (but not shape) to vary. 
A sky-map of the quantity 




A = 21n 



-£(S + B)- 
. ^(B) . 



(4) 



(the maximum likelihood ratio or MLR) is used to estimate and constrain the source location 
and the statistical significance of detection. The most probable source position corresponds to 
the maximum A = A in the sky-map, which formally obeys the probability distribution with 
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three degrees of freedom (source intensity and position). Thus, A > 14.2 indicates a statistically 

significant 3a detection. For such significant detections the source location is constrained by the 
quantity AA = A — A, which is distributed as X2 ^° that formal la, 2a and 3a statistical confidence 
regions in the sky-map arc defined where AA < 2.3, 6.2 and 11.8, respectively. 

The formal relation between likelihood ratio and is based on the implicit assumption of an 
accurate PSF and background model. The form of the background model is thus an important 
part of burst imaging, even though the number of background events during bursts is typically 
small. The dataspace distribution of background events present during GRBs is difficult to 
model due to complicated temporal and spatial variations. An empirical model provides the 
most practical means of approximating this, but background count rates are too low to provide 
a sufficient statistical sample of the full 3-d distribution. Fortunately, the (x,V') component of 
this distribution is reasonably approximated by the geometry function Q mentioned above. The 
ip component embodies much of the background complexity. It is determined independently for 
each burst by directly sampling events that were detected under orbital history conditions similar 
to the burst observation. In practice, this is done by sampling events within a window of satellite 
orbital parameters chosen to approximate the orbital history of CGRO at the time of a burst, and 
thereby the background environment during the burst. The distribution of directly sampled 
background events is combined with the geometry function (properly normalized) to form a 3-d 
background model that is then used in the imaging process. An example of real and modeled 
background is shown in Figure 3. 



2.4. Systematic Imaging Errors 

The approximations used in the imaging process are all potential sources of systematic errors 
in the derived burst source parameters (i.e., soTircc significance and localization). The major 
sources of error are: a) the use of a single PSF for all bursts, regardless of source location or 
energy spectrum; b) statistical fluctuations in the PSF; c) the inherent accuracy of Monte Carlo 
simulations used to compute the PSF; and d) the accuracy of the background model. The effects 
of each of these approximations have been studied with the use of simulations like those used 
to determine the PSF. Unfortunately, since all combinations of potential causes of error cannot 
possibly be simulated, this technique serves only as a predictive estimate of the true errors. 

Simulated sources with various, burst-like energy spectra, distributed throughout the FoV 
and mixed with samples of real background were subjected to the full imaging process with a 
variety of PSFs. The distribution of offsets between known source parameters and those derived 
through imaging analysis of the simulated bursts yields an estimate of systematic errors. We 
find that approximation (a) has the largest effect on source location accuracy, especially for 
??c ^ 30°, but that the error is never more than 0.75°. The other potential causes of location 
error contribute only minimally. Overall, simulated sources are correctly localized by the MLR 
imaging technique to within an error of 0.5° (~90% confidence). We also find that A (indicating 
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Fig. 3. — An example of the spatial distribution of COMPTEL background events during a 
gamma-ray burst. Observed (doia points) and modeled {solid lines) distributions are shown summed 
along each of the three dataspace coordinates. 

detection significance) is correctly determined to within an error of approximately ~30%. The 
largest contributing cause in this case is approximation (d), since the shape of the background 
in the observed COMPTEL dataspace is difficult to model accurately. The largest errors in A 
are for sources located near the earth's limb, where atmospheric gamma rays contribute most to 
the background and are not approximated well by our model. Bursts located near the earth's 
limb also result in source photons backscattered off the atmosphere, which confuse the true source 
signature. 
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3. Observations 

3.1. Burst Search Strategy 

To identify potential COMPTEL gamma-ray bursts, we rely on the BATSE transient event 
triggering system. This system automatically detects statistically significant count rate (50-300 
keV) increases above background in the BATSE large area detectors (Fishman et al. 1989). 
Triggers can occur from gamma-ray bursts, solar flares, magnetospheric activity. Earth occultation 
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of strong celestial sources and terrestrial atmospheric events. The COMPTEL burst search 
strategy is to examine data only from the times of BATSE triggers that are classified as GRBs 
by the BATSE analysis team (Mallozzi et al. 1993) and whose BATSE locations satisfy ??c < 65°. 
This strategy is effective in identifying COMPTEL GRB candidates. An exception occurs during 
BATSE data accumulation and readout periods following all triggers (Meegan et al. 1996). The 
BATSE trigger system is disabled during these accumulation periods (usually lasting either 4 min 
or 8 min) and insensitive to new, weaker events during data readout intervals (usually 90 min). 
These gaps in coverage must be taken into account when computing the exposure to bursts, as 



explained in §4.1 



In the period between 1991 April 19 and 1995 May 31 there were 3608 BATSE triggers, 
1297 of which have been classified as GRBs. Of these, 397 have BATSE locations with -d^ < 65°. 
COMPTEL data are available for only 302 of the 397 burst candidates due to gaps in telemetry 
caused by the failure of CGRO data recorders in 1992. Data around the time of each candidate 
trigger (±15 min) have been manually examined to identify potential burst intervals. Time 
intervals showing the most evidence for significant counts above background were subjected to 
the maximum likelihood imaging process. In cases where time intervals of significant emission 
were ambiguous, the approximate full duration of the burst as measured by BATSE was used for 
imaging. MLR maps were computed over a grid of sky coordinates with 1° x 1° spacing and a 
final burst detection threshold was applied to the images by requiring A > 20 — formally indicating 
a 3.8fT point source detection. The possibility of systematic error in the determination of A is 
accounted for by using a conservative detection threshold. This ensures that only significant 
excesses (>3(t) are accepted. We expect, at the most, one spurious detection in the search of 302 
candidates. 



3.2. COMPTEL Gamma-Ray Burst Location Catalog 

Over the four years of observations, we have accumulated a catalog of 29 gamma-ray bursts 
that satisfy the A > 20 detection requirement. These bursts are summarized in Table 1. Each 
detection is identified by the year (yy), month (mm) and day (dd) of the burst in the form GRB 
yymmdd. The next three columns give the number, Truncated Julian Day (TJD), and seconds 
of the day of the corresponding BATSE burst trigger. This is followed by the value of A and the 
most probable COMPTEL burst location expressed in Spacecraft (i^c^V'c), Celestial (a, (5; epoch 
J2000.0) and Galactic (/,6) coordinates. Preliminary locations of many of the bursts have been 
published earlier, as indicated. This catalog includes two bursts in addition to those reported 
in the preliminary analysis of Kippen et al. (1996b). One of these bursts (GRB 920627) was 
mistakenly neglected in the preliminary search and the other (GRB 950522) was beyond the 
time-interval of the earlier analysis. 

Statistical uncertainty regions around the most probable burst locations have been estimated 
by computing smoothed (cubic spline) contours of the MLR sky-maps at levels corresponding 
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COMPTEL Gamma-Ray Burst Location Catalog 
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■^This burst was detected by BATSE during the data readout period following a solar flare. 

''Angular radius of a circle having the same area as the irregularly shaped COMPTEL la confidence region. 

REFERENCES. — (1) Schonfelder et al. 1991; (2) Winkler et al. 1992; (3) Varendorff et al. 1992; (4) Winkler et al. 1992; (5) Connors et al. 
1993a; (6) Connors et al. 1993b; (7) Winkler et al. 1993; (8) Varendorff et al. 1993; (9) Kippen et al. 1993a; (10) Ryan, Kippen & Varendorff 
1993; (11) Bennett et al. 1993; (12) Ryan et al. 1994a; (13) Kippen et al. 1994a; (14) Kippen et al. 1994b; (15) Kippen et al. 1994c; (16) 
Ryan et al. 1994b; (17) Hanlon et al. 1994; (18) Kippen et al. 1995a; (19) Kippen et al. 1995b; (20) Greiner et al. 1995; (21) Winkler et al. 
1995; (22) Kippen et al. 1995c; (23) Kippen et al. 1995d; (24) Kippen et al. 1995e; (25) Kippen et al. 1996b. 



- 12 - 



to la, 2a and 2>a confidence. Celestial coordinate sky maps of these contours are shown in 
Figure 4. In general, the localizations are irregularly shaped regions, making it difficult to compare 
uncertainties between different bursts. A convenient approximate measure of the statistical 
location error is the radius of a circle having the same area as the region enclosed by the MLR 
contours. The value of this quantity ((Jeff) for the la confidence level is given in Table 1. The 
statistical location error is usually inversely related to the number of source events. Thus, the 
error radius of a high-fluence burst like GRB 940217 (cTcfr = 0.32°) is much smaller than that of 
a weaker detection like GRB 910627 (cJcfr = 1.63°). This behavior is complicated, however, by 
the character of the off-axis telescope response, which introduces additional uncertainty in the 
azimuthal ((/?c) direction for all bursts at large zenith angle (??c ^ 30°). The average statistical 
location uncertainty, as approximated by crgfj, for the full COMPTEL burst sample is 1.13°. A 
circle with this angular radius is included in each of the plots of Figure 4 for comparison. Note 
that systematic errors are not included in the MLR localizations. 

3.3. Measured Systematic Burst Location Error 

As indicated in the sky-maps of Figure 4, 26 of the 29 COMPTEL bursts have been 
independently localized through the method of arrival-time analysis using the CGRO and Ulysses 
spacecraft as two nodes of an Interplanetary Network (IPN) baseline (Hurley et al. 1992). This 
technique constrains the direction of a burst to a thin annulus on the sky, centered at coordinates 
(^arc) ^arc) with typical angular radius i?arc ~ 50° and width a^ < 0.1° (Hurley et al. 1996). 
The perpendicular offset p between the most probable COMPTEL burst location (Z, b) and the 
corresponding IPN annulus is given by 

p = COS"-^ [sin h sin 6arc + cos h cos 6arc cos (/ - Zarc)] " ^arc , (5) 

where positive and negative values indicate whether the COMPTEL location is inside or outside 
the annulus, respectively. This quantity provides a convenient measure of the total COMPTEL 
location error (statistical and systematic) in one dimension. We have used the CGKO / Ulysses 
IPN annuli of 26 bursts to directly estimate the average COMPTEL systematic burst location 
error (cf. Graziani & Lamb 1996). 

We began by simulating many COMPTEL "catalogs" of burst positions, with each position 
randomly blurred according to the statistical error distributions described by the MLR maps. 
An additional systematic blurring error was applied to the positions in the form of a symmetric 
Gaussian probability function with 68% confidence radius asys- For each random burst position, 
the quantity psira was determined from equation (|5|), assuming that the corresponding IPN annulus 
passes directly through the centroid of the input (for the simulation) burst position. Uncertainty 
in the IPN annulus width was incorporated by assuming a Gaussian error distribution with 
standard deviation cTr. Using the Kolmogorov-Smirnov test, we computed the maximum deviation 
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Fig. 4. — These plots show the locations of each of the COMPTEL bursts in Celestial coordinates 
(epoch J2000.0). Smoothed contours obtained through maximum likelihood imaging enclose la, 
2a and 3a statistical confidence regions. Where available, IPN timing annuli (Hurley et al. 1996) 
are overlaid, with the 3a error in the annulus width indicated by two concentric arcs. For scale 
reference, a circle corresponding to the average la area of all the bursts (angular radius 1.13°) is 
shown in each plot {thick line). 
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Fig. 4. — Continued. 
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Fig. 4. — Continued. 



of the cumulative distribution of observed p from that obtained with 2000 blurred catalogs 
(Psim)- By varying cTgys, and repeating the procedure, we determined how large fJsys must be in 
order to minimize the deviation between the observed and blurred distributions. We found that 
(Tsys = 0.6° ± 0.4° yields the best-fit within a 68% confidence interval. The observed and best-fit 
differential distributions of p are shown in Figure 5. This measurement is entirely consistent with 
the value of CTgys = 0.5° discussed above, and estimated by independent means. For all subsequent 
analysis we took Ggys = 0.5° and assumed that it was a constant for all bursts. When this error 
is added in quadrature to the statistical error radius (o"eff), the average COMPTEL location 
uncertainty increases to 1.25°. 




P (deg) 

Fig. 5. — Distribution of the angle p between most-likely COMPTEL burst positions and IPN 
annuli (at closest approach) for 26 bursts. The solid histogram shows the measured distribution, 
while the dashed curve is the best fit model incorporating statistical COMPTEL and IPN errors 
and a constant COMPTEL systematic location error of 0.5°. 



3.4. Combined COMPTEL/IPN Localizations 

The COMPTEL burst localizations are reduced over an order-of-magnitude in area by 
incorporating the constraints provided by IPN measurements. The COMPTEL localizations 
were used to limit the probable extent of the IPN annuli in the following manner. First, the 
probability distributions obtained from the observed MLR maps were convolved with a Gaussian 
systematic blurring function with crgys = 0.5°. This provided, for each burst, a measure of the 
total COMPTEL location error distribution about the most-likely position. The convolved maps 
were then evaluated at 0.1° intervals along the corresponding IPN annuli — yielding probability as 
a function of angular distance along the annuli. The most probable combined COMPTEL/IPN 
location of a burst corresponds to the point of maximum probability along the arc and confidence 
intervals were obtained from the integral probability distribution. Uncertainty in the width of 
the combined locations was assumed to be dominated by IPN measurement errors that are small 
compared to COMPTEL errors. 



TABLE 2 

COMPTEL/IPN Gamma-Ray Burst Locations'' 



Burst 02000 ^2000 I b Gs. (Jr Space- 
Namc (degrees) craft'' 



GRB 


910425 


89, 


,78 


-22 


.30 


227, 


.99 


-21, 


,03 


0, 


.56 


0, 


,053 


C,U 


GRB 


910503 


87. 


,43 


+38, 


.62 


172 


.13 


+5, 


,73 


0, 


,29 


0, 


,010 


C,U 


GRB 


910601 


309, 


,92 


+32 


.33 


73 


.82 


-5, 


,67 


0, 


.36 


0, 


,011 


c,u 


GRB 


910627 


198, 


,63 


-4, 


.35 


313 


.91 


+58, 


,04 


1, 


.43 


0, 


,031 


c,u 






198, 


,83 


-4, 


.52 


314 


.23 


+57, 


,84 


1, 


,35 


0, 


,036 


U,P 


GRB 


910709 




























GRB 


910814 


343, 


,78 


+29 


.31 


94 


.54 


-27, 


,09 


0, 


,35 


0, 


,018 


c,u 






343, 


,78 


+29 


.32 


94 


.54 


-27, 


,08 


0, 


,37 


0, 


,019 


U,P 


GRB 


911118 


166, 


,16 


-23 


.21 


273 


.30 


+33, 


,46 


1, 


,98 


0, 


,017 


c,u 






166, 


,48 


-23 


.23 


273 


.62 


+33, 


,59 


1, 


,78 


0, 


,019 


U,P 






166, 


,06 


-22 


.59 


272, 


.84 


+33, 


,96 


0, 


,61 


0, 


,086 


-1- ^ 


GRB 


920622 


161, 


,42 


+46, 


.86 


165, 


.83 


+58, 


,44 


0, 


,35 


0, 


,009 


c,u 


GRB 


920627 


169, 


,04 


-8, 


.36 


266 


.57 


+47, 


,69 


1, 


,14 


0, 


,049 


c,u 


GRB 


920830 


266, 


,05 


-74, 


.67 


318 


.75 


-21, 


,82 


1, 


,08 


0, 


,018 


c,u 


GRB 


930118 




























GRB 


930131 


183, 


,65 


-9 


.72 


288, 


.10 


+52, 


,10 


0, 


,97 


0, 


,003 


c,u 


GRB 


930309 


324, 


,93 


+50, 


.72 


94 


.90 


-1, 


,43 


0, 


,66 


0, 


,035 


c,u 


GRB 


930612 


106, 


,88 


-70 


.88 


281 


.94 


-24, 


,21 


0, 


,98 


0, 


,053 


c,u 






110, 


,65 


-69 


.97 


281, 


.34 


-22, 


,75 


1, 


,05 


1, 


,289 


C,M 


GRB 


930704 


99, 


,46 


+62, 


.98 


152 


.32 


+22, 


,48 


1, 


,53 


0, 


,113 


c,u 






99, 


,99 


+62 


.40 


153 


.00 


+22, 


,55 


1, 


,05 


2, 


,319 


C,M 


GRB 


931229 




























GRB 


940217 


29, 


,50 


+3, 


.83 


152 


.97 


-55, 


,11 


0, 


,29 


0, 


,020 


c,u 


GRB 


940301 


102, 


,01 


+64, 


.20 


151 


.35 


+23, 


,88 


0, 


,37 


0, 


,090 


c,u 


GRB 


940314 


355, 


,74 


-2 


.56 


86 


.35 


-60, 


,49 


1, 


,29 


0, 


,220 


c,u 


GRB 


940520 


324, 


,71 


+7, 


.39 


62 


.34 


32, 


,05 


0, 


,66 


0, 


,059 


c,u 


GRB 


940619 


298, 


,59 


-30 


.05 


10 


.81 


-25, 


,95 


0, 


,91 


0, 


,081 


c,u 


GRB 


940708 


300, 


,59 


+23 


.27 


61 


.65 


-3, 


,98 


0, 


,64 


0, 


,038 


c,u 


GRB 


940728 


85, 


,94 


-38 


.38 


243 


.80 


-29, 


,00 


0, 


,74 


0, 


,131 


c,u 


GRB 


940921 


207, 


,01 


-17, 


.24 


321, 


.77 


+43, 


,56 


0, 


,55 


0, 


,121 


c,u 


GRB 


941017 


298, 


,14 


+13, 


.09 


51, 


.68 


-7, 


,22 


3, 


,16 


0, 


,016 


c,u 


GRB 


950208 


344, 


,78 


+55, 


.94 


107 


.63 


-3, 


,57 


2, 


,95 


0, 


,022 


c,u 


GRB 


950421 


69, 


,78 


-61, 


.77 


272, 


.19 


-39, 


,14 


1, 


,14 


0, 


,226 


c,u 


GRB 


950425 


164, 


,69 


-33 


.77 


277, 


.62 


+23, 


,52 


0, 


.66 


0, 


,207 


c,u 


GRB 


950522 


111, 


,82 


+19 


.22 


199 


.32 


+16, 


,35 


0, 


.57 


0, 


,063 


c,u 



''IPN data from Hurley et al. (1996). 

''IPN spacecraft: C=Compton(BATSE), \]=Ulys8es, F=Pioneer Venus 
Orbiter, M=Mars Observer. 
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Combined COMPTEL/IPN burst localizations are listed in Table 2, where the most probable 
combined locations are expressed in Celestial (epoch J2000.0) and Galactic coordinates. The next 
column gives the mean 68% confidence error cXa in the angular distance along the arc, as measured 
from the most probable position. This is followed by the uncertainty in the arc width Ur, as 
determined from IPN timing measurements (Hurley et al. 1996). Those bursts that were observed 
by more than two spacecraft are indicated. The average COMPTEL/IPN burst localization {la 
confidence) corresponds to a short, thin arc-segment with total length 2 x (a"a) = 1.97° and total 
width 2 X (cji.) = 0.13°. For one burst (GRB 911118), IPN measurements from three spacecraft 
provide a nearly unambiguous determination of the source location where the individual annuli 
intersect. The IPN location of this event is within the COMPTEL statistical la MLR contour 
(see Figure 4). 

4. Spatial Distribution Analysis 

Despite their limited number, the COMPTEL GRB catalog locations provide a useful measure 
of the global angular distribution of GRBs. In principle, the catalog can also be used to probe 
the GRB distance scale through the distribution of burst intensities (e.g., Pendleton et al. 1996; 
Fcnimore et al. 1993). However, the COMPTEL measurements are ill-suited for this purpose 
because of severe observational biases — most notably the strong dependence of burst sensitivity 
on fiuence (and duration), the large variation of sensitivity across the FoV and the inability to 
measure intense burst episodes due to deadtime. These biases would require large, uncertain 
correction factors at nearly all observed intensity levels. The COMPTEL locations can, however, 
provide some insight into the burst distance dilemma through correlations with known objects 
and by providing a lower limit on the distance to individual events through comparison with IPN 
measurements. 

We analyzed the distribution of COMPTEL burst locations using two inclusive samples 
chosen for easy comparison with the results of earlier analyses. The three-year sample includes 18 
bursts observed between 1991 April 19 and 1994 March 1, while the four-year sample incorporates 
the full catalog of 29 bursts. Galactic coordinate distributions of the COMPTEL statistical 
location constraints for both burst samples are shown in Figure 6. The null hypothesis is that 
GRB sources have an isotropic angular distribution on all size scales, with no preferred direction 
or clustering. We use the COMPTEL bursts as a sample of the full GRB population to test this 
hypothesis. The small number of bursts means that statistical fluctuations are the dominant 
source of error and that this error is non-Gaussian. Statistical tests are therefore conducted by 
Monte Carlo methods. The major observational biases present are the non-ideal burst locations 
and uneven sampling of the sky, both of which affect the interpretation of the data. 
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+90 +90 




Fig. 6. — The COMPTEL statistical {la, 2a and 3a) location contours (top) and corresponding 
burst exposure functions (bottom) plotted in Galactic coordinates. The leftmost plots correspond 
to the three- year observation period (18 bursts), while the plots at right include the full four-year 
interval (29 bursts). 



4.1. Burst Exposure Function 

The limited extent and varying sensitivity of the COMPTEL FoV, combined with CGRO's 
low Earth orbit and sporadic pointing schedule result in a complicated sampling pattern over the 
sky that modifies the intrinsic angular source distribution. We must determine this pattern — the 
burst exposure function — in order to correct the observational bias it produces. 

Assuming that bursts occur isotropically in space and uniformly in time, the relative 
probability Pexp of detecting a burst from a particular direction in the sky compared to any other 
direction is given by 



Pexp — ^on 



Js^i^ dS 



where ton is the time spent observing the direction of interest, Smin is the minimum fluence for a 
detection and ^ is the differential burst size-frequency distribution (i.e., the number of bursts 
per year with 0.75-30 MeV fluence between 5 and S + dS). Evaluation of Pgxp is complicated 
by the function S'minj that depends on the source position within the FoV, the background rate 
and the burst energy spectrum. We approximated Pexp by ignoring background and taking the 
burst energy spectrum to be the average power law (^-2-6; Kippen et al. 1996a, 1997)— leaving 
5inin as a function only of sky-position through equation (|2|). A further problem is that is 
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not well-known in the COMPTEL energy range. For simplicity, we used a power law of the 
form ^ cic S~^/^. For strong bursts, such as those observed by COMPTEL, this is a reasonable 
approximation (Fenimore et al. 1993). 

The above approximations were used to compute Pexp for every 16 seconds of observation on 
a grid of sky coordinates with 2° x 2° spacing. Directions occulted by the earth were excluded, as 
were all epochs when COMPTEL and/or the BATSE trigger system were disabled or the spacecraft 
was not telemetering data. Directions outside the -de < 65° search window were also excluded. 
Sky-maps of the burst exposure function accumulated over the two analysis intervals are shown 
in Figure 6. The maximum value of Pexp in the three-year (four-yeax) exposure map is 426 (624) 
effective days. In both maps, the minimum exposure is at least 30% of the maximum — meaning 
that all portions of sky have been observed at a reasonable level. The smooth fluctuations on top 
of this minimum exposure reflect the frequent Earth occultation of directions along the Celestial 
equator and the schedule of the CGRO pointing program. In particular, the south Galactic polar 
region received the least exposure, while a great deal of time was spent observing the Galactic 
plane and Virgo (Z = 270°, b = 60°) regions. The effects of different spectral shapes, burst intensity 
function and the inclusion or exclusion of BATSE data accumulation periods are small (<10%) 
compared to the overall magnitude of the structures in the exposure maps. 



4.2. Large-Scale Angular Distribution 

Studies of large-scale angular structure have been most concerned with the question of 
whether burst sources reside inside or outside our Galaxy. Two convenient and sensitive statistical 
tools for this particular purpose are the dipole and quadrupole moments in Galactic coordinates 
(e.g., Hartmann & Epstein 1989; Paczyhski 1990; Briggs 1993). The dipole moment (expressed 
as (cos O) , where O is the angle between a burst and the Galactic center) measures the amount 
of concentration towards the center of the Galaxy, while the quadrupole moment (expressed as 
(sin^ 6 — |), where b is the Galactic latitude of a burst) indicates the amount of concentration 
towards the Galactic plane. Significant non-zero values of the moments indicate deviations from 
an isotropic source distribution. Similar statistics, expressed in Celestial coordinates (i.e., (cos0), 
where B is the angle from the Celestial origin and (sin^ (5 — |), where S is declination), are more 
sensitive to the biases caused by the non-uniform COMPTEL exposure function. 

The dipole and quadrupole moments expected for an isotropic angular distribution of sources 
are modified by the non-uniform burst exposure function. The magnitude of this effect was 
estimated by simulating random catalogs of burst locations chosen according to the COMPTEL 
exposure function maps. An ensemble of moments computed from 5 x 10'^ simulated catalogs 
provides a measure of the exposure-induced deviation from isotropy and the statistical uncertainty 
inherent in the number of bursts in each catalog. The largest deviation from the isotropic value 
appears in the Celestial quadrupole moment, where the exposure function produces a 0.7a bias 
towards polar declinations (both three and four-year samples). The exposure-induced bias in 
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Galactic coordinates is much smaller (see examples shown in Figure 7). 




-\ — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — p 




(cos 0> 

Fig. 7. — Simulations of Celestial quadrupole moment (a) and Galactic dipole moment (b) 
distributions for 29 burst locations. The narrowest distributions in each case represent the values 
observed by COMPTEL, including the effect of location errors. The remaining curves show the 
statistical distributions expected from isotropy assuming uniform {dashed curves) and non- uniform 
{solid curves) burst exposure functions. 

Observed moments were estimated by randomly blurring each of the measured burst locations 
according to the total (statistical and systematic) location error distributions. An ensemble 
of moments computed from many blurred catalogs provides a measure of the error induced by 
imprecise locations. As illustrated in Figure 7, this effect is small compared to the statistical 
uncertainty caused by the low number of bursts. The difference between the observed and 
expected moments yields the true or corrected deviation from isotropy. Corrected values of the 
dipole and quadrupole moments are given in Table 3 in both Galactic and Celestial coordinates. 
Uncertainties in the corrected moments reflect 68% confidence limits on the difference between 
the data and the simulations. Within the sizable uncertainties, there are no significant deviations 
from large-scale isotropy evident in the COMPTEL burst samples in either Galactic or Celestial 
coordinates. 
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TABLE 3 

Large-Scale Angular Distribution Measures 



Statistical 
Test 


Three- Year Sample 


Four- Year Sample 


Value 


Deviation^ 


Value 


Deviation^ 


Galactic 










(cos e) 


-0.124 ± 0.140 


-0.886 


-0.016 ± 0.110 


-0.145 




+0.001 ± 0.068 


+0.015 


-0.023 ± 0.053 


-0.434 


Celestial 










(cos e) 


-0.091 ± 0.129 


-0.705 


+0.026 ± 0.101 


+0.257 


(sin2 S-D 


+0.034 ± 0.073 


+0.466 


-0.022 ± 0.057 


-0.386 



^Number of standard deviations of the measured value from isotropy. 

4.3. Small-Scale Angular Distribution 

The dipole and quadrupole moments are sensitive only on the largest of angular scales 
{9 > 60°). Different statistical tests are necessary to investigate the small-scale angular distribution 
of burst sources. The two standard approaches are the angular autocorrelation function w{9) and 
the nearest neighbor function nn{9). Both statistics are sensitive to small-scale angular structure 
on the order of the instrumental resolution, which could suggest either intrinsic source clustering or 
repeated bursts from individual sources. It has been shown that the angular correlation function 
provides superior sensitivity to small-scale angular structure when the number of sources is large, 
such that the average separation between neighboring positions is less than the average location 
accuracy (Brainerd 1996). This is not the case for the small sample of well- localized COMPTEL 
bursts, where the average neighbor separation (~20°) is much greater than the average location 
accuracy (~1.25°). Hence, both tests were used. 

The angular correlation function, first applied to GRB positions by Hartmann &: Blumenthal 
(1989) and Hartmann, Linder &: Blumenthal (1991), measures the probability of finding a pair of 
sources separated by an angle 9. It is computed by comparing the number of source pairs A^obs) 
with angular separations in the range {9,9 + d9), to the number of pairs iVexp expected from a 
random sample of an isotropic distribution with 

wi9) = ^-l. (7) 

'exp 

The presence of small-scale angular clustering is evident in w{9) as a significant positive deviation 
from zero at small angles. The nearest neighbor function im{9), first used on GRB positions by 
Quashnock k. Lamb (1993), differs slightly in that it compares the number of sources with one or 
more neighbors separated by {9, 9 + d^) to the number expected in a random, isotropic sample. 
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Small-scale clustering is evident in nn{6) as a significant excess above the expected value at small 
angles. Both statistics are affected by imprecise source locations, which spread the intrinsic signal 
of clustering over a range of 9 characteristic of the average location accuracy. The signal of 
small-scale clustering is maximized by choosing bins in 6 with the same characteristic size. 

Figure 8 shows the angular correlation and nearest neighbor distributions for the two samples 
of COMPTEL burst locations. Angular bins of 2° were used to approximate the average location 
accuracy at the 90% confidence level, where the signal from clustering would be most apparent. 
The distributions of A^exp expected from isotropy were determined by taking the average of 
5 X 10^ simulated catalogs, each chosen at random according to the corresponding burst exposure 
function. The effect of non-uniform exposure on the A'exp distributions is small compared to 
statistical fluctuations among the random catalogs (< O.lcj). The statistical significance of 
observed deviations in excess of the isotropic expectation was estimated by the fraction of 
simulated catalogs Q with an equal or greater excess in the angular bin in question. 
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Fig. 8. — Angular auto-corrclation functions w{0) and nearest neighbor distributions rn{9) of 
two samples of COMPTEL burst locations. Each plot compares the distribution measured with the 
most-likely locations {solid histograms) to the distributions expected from the non-uniform sampling 
of isotropy (solid curves). Dashed curves indicate the magnitude of statistical deviations from the 
expected distributions at the 90% significance level. 
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The most significant small-angle excess is evident in the first 2° bin in the three-year 
catalog distributions, where Q = 4.8 x 10~^. This excess occurs because two of the bursts 
(GRB 930704 and GRB 940301) are localized to the same direction within their combined location 
errors — something that is not expected in a random sample of only 18 sources (Ryan et al. 
1994a). This near coincidence of two bursts occurring eight months apart from one another was 
investigated by Kippen et al. (1995b) and Hanlon et al. (1995), who computed a similar chance 
probability for the observation {Q = 3.0 x 10~^). The value of Q computed here is somewhat 
larger due mainly to the addition of the previously neglected burst GRB 920627. Kippen et 
al. (1995b) could neither prove, nor rule out the possibility that the two bursts originated from 
a single source — even when all other available measurements were included (IPN, BATSE and 
CGRO-EGRET) . In the larger four-year sample the significance of the excess at ^ < 2° is reduced 
to Q = 1.2 X 10~^ — a result of the fact that no additional bursts from the possible repeating 
source were observed in the year following GRB 940301, and that no other coincident pairs were 
measured elsewhere. We therefore conclude that the four-year sample is consistent with isotropy 
on small scales. 



4.4. Correlation with Galaxy Catalogs 

Examination of correlations between GRBs and catalogs of known objects is similar to 
the study of self-clustering described above. Statistical tools like the nearest neighbor test and 
the angular correlation function are simply modified by computing angles between GRB/object 
pairs, rather than angles between pairs of different bursts. As in the case of self-clustering, cross 
correlations between GRBs and known objects are affected by burst location errors, which spread 
the cross-correlation signal over an angular range characteristic of the average burst location 
error. Although many different types of astronomical catalogs have been examined, we chose to 
concentrate only on galaxy clusters and radio-quiet quasars, which are suggested to correlate with 
BATSE SB burst locations (Kolatt & Piran 1996; Marani et al. 1997; Schartel, Andernach & 
Greiner 1997). In particular, we used the galaxy cluster catalog of Abell, Corwin k, Olowin (1989, 
hereafter AGO), and the quasar catalog of Veron-Getty & Veron (1996, hereafter VGV). The full 
AGO catalog (AGO-1) contains 5250 galaxy clusters, categorized by several parameters, including 
richness class R and distance class D. A subset of the 185 AGO galaxy clusters with R>1 and 
D < 4 (AGO-2) was examined separately since this particular subset was suggested to show a 
correlation with the most accurate BATSE burst locations (Marani et al. 1997 and references 
therein). The VGV catalog contains 11,662 quasars and active galactic nuclei. We examined only 
the radio-quiet quasars, as selected by Schartel, Andernach & Greiner (1997). The full radio-quiet 
quasar sample (VGV-1) contains 7146 objects. A subset of 967 radio-quiet quasars with redshift 
z < 1.0 and absolute magnitude Mabs < —24.2 (VGV-2) was also examined since this sample 
shows the strongest correlation with BATSE GRBs. 

The computation of angular cross-correlation between GOMPTEL GRBs and AGO/ VGV 
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Fig. 9. — Angular cross-correlation function w{d) between 29 COMPTEL burst locations and: (a) 
185 AGO galaxy clusters with richness class R > 1 and distance class D < 4; (b) 967 VCV radio- 
quiet quasars with rcdshift z < 1.0 and absolute magnitude Afabs < —24.2. Each plot compares 
the distribution observed with the most-likely COMPTEL burst locations {solid histograms) to the 
distribution expected, on average, from random samples of the non-uniform COMPTEL exposure 
function (solid lines). Dashed curves indicate the magnitude of statistical deviations from the 
expected distributions at the 90% significance level. 





TABLE 4 

GRB/Galaxy Correlation Statistics 



Catalog 


Objects 


^obs 


^exp 


Q 


ACO-1 


5250 


46 


42.9 ± 12.0 


0.390 


ACO-2 


185 


3 


1.6 ± 1.4 


0.218 


VCV-1 


7146 


25 


53.8 ± 45.0 


0.736 


VCV-2 


967 


10 


7.7 ± 5.7 


0.286 
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catalog locations was performed as in the case of self-clustering. In each 2° angular bin the number 
of observed GRB/catalog pairs A'obs was compared to the average number obtained from 5 x 10'* 
random COMPTEL "catalogs" A'exp (sampled from the four-year COMPTEL burst exposure 
function) through equation (^). Table 4 lists the results for GRB/catalog pair angular separations 
9 < 2° . The uncertainty in A'^exp given in the table reflects the standard deviation of the random 
catalogs and Q is the statistical significance of the observed number of pairs compared to the 
expected number. None of the four catalog samples show a significant correlation with COMPTEL 
GRB locations at small angles, where one expects the signal if an association exists. Angular 
cross-correlation functions for the ACO-2 and VCV-2 object catalogs are shown in Figure 9. 



4.5. COMPTEL/IPN Burst Distance Limits 

While a "local" GRB source population has not been examined in as much detail as the 
more distant possibilities, such an origin (i.e., sources within 10^ AU) is not inconsistent with 
the angular and intensity distributions observed by BATSE (Maoz 1993; Horack et al. 1994). 
Lower limits on the distance to burst sources are therefore important. As pointed out by Hurley 
(1982), timing differences between four or more widely separated IPN spacecraft in a favorable 
configuration could be used to yield a lower limit on the burst source distance. While only a few 
events have been observed by this many IPN spacecraft, the idea can also be applied to the bursts 
localized by COMPTEL and the two-spacecraft, CGRO/ Ulysses IPN. 

If a burst is observed by two spacecraft separated by a distance 6, the angle a between the 
burst direction and the inter-spacecraft baseline is 



where c is the speed of light, D is the distance to the burst source and At is the observed difference 
in burst arrival times between the two spacecraft (Connors et al. 1993a). Thus, cos a is larger 
than that inferred from the time-delay alone (cos a = cos -Rare = cAt/b) by a parallax factor that 
is inversely related to the source distance. If the offset of a burst position from an IPN annulus 
were measured, it would reveal the source distance. The fact that COMPTEL burst localizations 
agree with IPN annuli to within the COMPTEL location accuracy places a lower limit on the 
distance to burst sources. The severity of this limit improves with the accuracy of the COMPTEL 
localization (in the direction perpendicular to the IPN annulus) and the length of the IPN baseline 
(b). 



According to the procedure described in §3.^^, the total (statistical and systematic) COMPTEL 



localization of a burst can be maximized along an annulus to reveal the most-likely combined 
position and the probability at that position. In computing a quantitative COMPTEL/IPN 
distance limit this technique was applied to a set of concentric annuli with angular radii a 
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TABLE 5 

COMPTEL/IPN Gamma- Ray Burst Distance 
Lower Limits^ 



Burst 


-^arc 


CTr 


b 




Name 


n 


n 


(AU) 


(AU) 


GRB 910601 


52.936 


0.011 


3.708 


109.4 


GRB 920622 


41.335 


0.009 


5.800 


79.8 


GRB 920830 


86.268 


0.018 


6.271 


180.3 


GRB 940217 


69.334 


0.050 


3.038 


137.5 


GRB 940921 


54.361 


0.121 


2.583 


90.1 



IPN data from Hurley et al. (1996). 



progressively smaller than Rave- This yielded — with the use of equation (^) — the COMPTEL 
location probability as a function of source distance D. Uncertainty in the width of the IPN 
annuli was incorporated by substituting R^rc + Cr in the calculation of probability as a function 
of D. Integration of this function at the desired confidence level yields a lower limit. The most 
constraining distance lower limits (-Dmin) at the 2a confidence level are listed in Table 5. Even 
when all known systematic location errors are conservatively included, the distance to these bursts 
must be at least ~100 AU in order to be consistent with the COMPTEL/IPN measurements. 

5. Summary and Discussion 

We have demonstrated the ability of COMPTEL to accurately localize GRB sources through 
direct Compton scatter imaging. This ability has been used to accumulate a unique collection 
of locations for bursts having a high fluence of gamma rays in the energy range 0.75-30 MeV. 
Analysis of the COMPTEL locations has shown that these bursts, apparently like all others, are 
consistent with an isotropic angular distribution of sources. Furthermore, we find no statistically 
significant excess of GRB locations clustered at small angles, which would be expected if many 
of the sources we observed burst repeatedly on time-scales of several months. We also find no 
significant cross-correlation with AGO galaxy clusters or radio-quiet VCV quasars, and have 
placed a lower limit on the distance to GRB sources by comparing the COMPTEL locations to 
IPN timing annuli. The COMPTEL measurements thus contribute to our understanding of the 
GRB phenomenon in each of the three currently feasible source distance regimes. 

For a local heliospheric origin, the distance constraints required by COMPTEL/IPN burst 
location measurements dispute any scenarios in which the sources of brighter (and therefore closer, 
assuming bursts are roughly "standard luminosity candles") GRBs are nearer than ~100 AU 
(approximately the distance of the solar wind termination shock) . Relativistic dust grain models 
(Grindlay & Fazio 1974; Dasgupta 1979) and nearby (40-400 AU) colliding comet scenarios (White 
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1993, 1994) are therefore unfavorable. More distant Oort cloud models involving mergers between 
comets, antimatter and/or primordial black holes (Katz 1993; Dermer 1996; Bickert & Greiner 
1993) are not excluded by our distance limits. However, dynamical simulations of potential Oort 
cloud GRB progenitor distributions indicate that it is difficult to restrict GRBs only to the outer, 
more distant portions of the cloud (Clarke, Blaes k. Tremaine 1994). In order to obey the angular 
and intensity distributions measured with BATSE and the distance lower limits reported here, 
Oort cloud models must invoke an ad hoc assumption that the GRB progenitor distance has a 
sharp cutoff between 100-500 AU (e.g., Horack et al. 1994). Future instruments capable of more 
accurate GRB localization (e.g., BcppoSAX and HETE-2) may increase the distance limit. For 
example, a burst localized with arc-minute accuracy could (when compared to an IPN annulus) 
raise the distance lower limit to thousands of AU. 

The Galactic dipole and quadrupole moments of COMPTEL burst locations independently 
confirm the long-known result that strong bursts are distributed isotropically (Atteia et al. 1987; 
Hartmann k, Epstein 1989). In addition, the COMPTEL measurements mean that GRBs with 
predominantly hard spectra are probably from the same angular source distribution as the main 
bulk of events observed by BATSE, even though the radial distribution of soft and hard events 
measured with BATSE intensities may differ (e.g.. Belli 1995; Kouveliotou et al. 1996). These 
results are important in the context of Galactic halo scenarios, where one expects all the nearer, 
brighter bursts to exhibit a significant positive dipole moment due to the 8.5 kpc offset of the Sun 
from the Galactic center. Unfortunately, the absence of a dipole in the COMPTEL bursts is not 
constraining due to the large uncertainties involved. 

In the cosmological distance regime, the absence of significant small-angle correlations among 
the COMPTEL locations is important, since it bears upon the question of burst repetition. Given 
that the larger four-year sample of bursts does not show any significant evidence for small-scale 
clustering, we must conclude that the two-burst overlap between GRB 930704 and GRB 940301 
was either a statistical fluke, or that the two bursts were from a rare repeating source (or source 
cluster), whose rate of repetition at the COMPTEL sensitivity level is low. If the latter case is 
true, the lack of excess BATSE bursts from the GRB 930704/GRB 940301 region indicates that 
the repetition rate from this source is also low at the BATSE sensitivity level. The rarity of 
repeating GRB sources can be estimated from the COMPTEL observations, where the fraction 
of observed repeating burst sources must be less than 2/29 = 7%. We can obtain this number 
without the use of complicated statistical tests, as are required in the analysis of BATSE locations 
(e.g., see Meegan et al. 1995; Tegmark et al. 1996), since any other possible repeater candidates 
would be immediately obvious given the small number of COMPTEL bursts and the accuracy of 
the locations. It is interesting that limits on the fraction of repeating burst sources obtained with 
BATSE (incorporating thousands of bursts) are about the same as the 7% given here, even though 
BATSE detects bursts that are an order of magnitude weaker (Meegan et al. 1996; Tegmark et al. 
1996). It thus appears that either: a) repeating burst sources are rare, b) the typical repetition 
time-scale is >1 yr, or c) repeated bursts are too weak to be reliably detected, even by BATSE. 
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More constraining results will come from continued observations and from future instruments with 
improved location accuracy. 

Of further importance in the cosmological distance regime is the lack of significant correlation 
between COMPTEL bursts and cither AGO galaxy clusters or radio-quiet quasars. These findings 
question the reality of BATSE/ACO/VCV associations that have been suggested. According to 
Marani et al. (1997) the strongest cross-correlation signal for BATSE/ACO clusters (~3.5-4a 
significance) is obtained for the 27 BATSE 3B bursts with the smallest location errors and that 
the signal weakens as bursts with larger location errors are included. This could indicate that 
bright (or high fluence) bursts are more highly correlated with AGO galaxy clusters than weak 
events, since bright bursts tend to yield smaller statistical location errors than weaker bursts. This 
behavior is exactly what we would expect of cosmological GRBs, where only the closest, brightest 
events would be correlated with relatively nearby AGO galaxies. However, the same correlation 
should exist in the GOMPTEL bursts, since they are the same average brightness and fluence as 
the 27 well-localized BATSE 3B events. In fact, the two samples share six bursts in common. 
The only signiflcant difference between the GOMPTEL and BATSE burst samples is that the 
GOMPTEL events have harder spectra, on average, than the BATSE bursts. It is possible that 
hard events are not as strongly correlated with AGO galaxy clusters as arc soft bursts, but it is 
more likely that the BATSE/AGO correlation is a statistical fluke. This conclusion is supported by 
Hurley et al. (1997), who flnd that a majority of the AGO clusters that are correlated with BATSE 
locations are inconsistent with IPN timing annuli — implying that the BATSE/AGO correlation 
is statistical in origin. The lack of a correlation between GOMPTEL burst locations and VGV 
radio-quiet quasars is less meaningful, since the number of GOMPTEL bursts is much smaller than 
the number of BATSE events (80 bursts with location errors < 1.7°) used by Schartel, Andernach 
&; Greiner (1997). The weak GRB/quasar correlation using the more numerous BATSE locations 
may be below the sensitivity of the GOMPTEL sample due to its larger statistical uncertainty. 

While the collective analysis of GOMPTEL burst locations offers some insight into the GRB 
mystery, perhaps the greatest value of the locations is in facilitating rapid counterpart search 
efforts. The recent discovery of x-ray and optical emission fading on a time-scale of hours/days 
from the BeppoSAX locations of GRB 970228 and GRB 970508 (Gosta et al. 1997a, 1997b; van 
Paradijs et al. 1997; Bond 1997) highlights the value of rapid follow-up studies using accurate 
burst locations. During the course of the GGRO mission, we have developed a rapid response 
system to quickly compute GOMPTEL GRB locations and distribute them to multi-wavelength 
observers (primarily optical and radio) around the world (Kippen et al. 1993). The advantage 
of using the relatively accurate GOMPTEL locations is that they allow deeper searches than 
previous, and most other current efforts. The system has been used for several bursts, resulting 
in some of the most sensitive radio/optical observations of burst locations ever obtained within 
several hours (see e.g., Harrison et al. 1995; McNamara et al. 1996 and references therein). While 
no fading counterparts have been detected as a result of the GOMPTEL rapid response effort, the 
observations obtained point the way for future studies, which must concentrate on deep exposures 
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within minutes of a GRB. This is the goal of continuing COMPTEL rapid response efforts, as well 
as several on-going and planned new missions. If the optical lightcurves of GRB 970228 and GRB 
970508 are typical then wide-field instruments guided by a COMPTEL GRB location within ~1 
hour may detect the fading counterpart. 
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